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The baryonic acoustic oscillations are features in the spatial distribution of the galaxies which, 
if observed at different epochs, probe the nature of the dark energy. In order to be able to 
measure the parameters of the dark energy equation of state to high precision, a huge sample 
of galaxies has to be used. The Large Synoptic Survey Telescope will survey the optical sky 
with 6 filters from 300 nm and 1100 nm, such that a catalog of galaxies with photometric 
redshifts will be available for dark energy studies. In this article, we will give a rough estimate 
of the impact of the photometric redshift uncertainties on the computation of the dark energy 
parameter through the reconstruction of the BAO scale from a simulated photometric catalog. 

1 Introduction 

Before the time of recombination, the competition between the photon pressure that tends to 
drive the baryons away from the dark matter potential wells, and the gravitational attraction 
between the baryons and the dark matter, induced acoustic waves in the photon-baryon plasma. 
We call these acoustic waves the Baryonic Acoustic Oscillations (BAO). When the photons 
and the baryons are decoupled the oscillations forze until the baryon-drag epoch. The CMB 
temperature anisotropy peaks are directly related to the sound horizon (r s (zd) = 100h _1 Mpc 
comoving) at the decoupling epoch (z^ = 1100). Similarly, the correlation function of the 
galaxies, that is supposed to trace the matter density field, exhibits a bump at the BAO scale, 
which, by taking the Fourier transform, translates into oscillations (kBAo) i n the matter power 
spectrum. As time evolves, the Universe expands according to the equations of state of its 
different constituents so the angular diameter distance as a function of redshift, which is related 
to the comoving BAO wavelength, also evolves. We propose a simulation that should allow an 
estimation of the uncertainties of the dark energy equation of state parameters as they should be 
measured with the telescope LSST^, by reconstructing Ubao a t different redshifts and estimating 
its uncertainties. 



2 Simulation pipeline 



2.1 Fiducial cosmology 

We start by setting the value of the cosmological parameters, roughly the same as the WMAP-5 
values for a ACDM cosmology: the dark matter density fi m = 0.27, the dark energy density 
r^A = 0.73, the baryon density = 0.04, the Hubble parameter h = 0.71, the dark energy 
equation of state w = —1, the spectral index n = 1 and the normalisation of the density 
fluctuations = 0.8. 

2.2 Matter power spectrum, density field and galaxy number 

The matter power spectrum for a redshift z, P(k, z) is computed with a model linear theory for 
the fiducial parameter values, from the fitting formulae of Eisenstein et al^ and is scaled to z 
using the approximation for the growth function G(z) from Carroll et al^ 

#(a) f a da' 
P{k,z) = G{z) 2 P(k,0), 

where P{k, 0) is the matter power spectrum measured today. A cuboid of Fourier coefficients 
is generated from PQs, z) then Fast-Fourier Transformed into a cuboid of over-densities, S(r) = 
(p(r) — p)/p, where p(r) is the matter density, and p is the background mean matter density. 
The matter density field relates to the number of galaxies in a cell of center r and volume V c 
with: iVj(r) = (<5(r) + 1)V C x n where n is the mean number density of galaxies per unit volume. 
We determine this density n by integrating the luminosity function (LF) (j>(M) that gives the 
probability of the expected number density of galaxies per unit volume and unit of absolute 
magnitude M: 

n = J (p(M)dM. 

The LF is a Schechter parametric function, and the values of the parameters are the ones 
measured by the GOODS survey from Dahlen et aPl It is, in fact, a function of the type of the 
galaxy: early, late and starburst and is measured in the B filter of the GOODS survey. 

2.3 Observed galaxy catalog 

Each of the simulated galaxies is characterized by an absolute magnitude Mb and a spectral 
energy distribution (SED), that is fixed by the LF. The SED is taken to be one of 50 SED that 
are interpolated among the four CWW templates from Coleman et al^ and two starburst SED 
from Kinney et al^H The apparent magnitudes in the 6 LSST bands (u, g,r,i, z,y) are then 
computed via: 

m x = M B + p(z) + Kbx(z), 

where X is a LSST filter, p is the distance modulus, z the galaxy redshift, and Kbx is the K- 
correction that adjusts the magnitude from the rest-frame B band of GOODS to the observed- 
frame X band. The K-corrections are computed following the formulae in Hogg et al^. A 
Gaussian photometric uncertainty is added to mx and is equal to: 

o-(m x ) = \J (0.04 - 7 )io°- 4 (™x-m 5 ,x) + 7 io0.8(m Jf -m 5 ,x) ) 

where m§x is the 5a detection magnitude in the X band. This equation and the form of m^x 
are given in the LSST Science BookH. A systematic uncertainty is added in quadrature, and is 
set to 0.005. If mx > m^^x the galaxy is not observed in this band. 



2-4 Photometric redshift 



For each galaxy, the photometric redshift is computed using our own template fitting method 
based on a chi square minimization, that is very similar to the public code LEPHARE or 
HYPERZE. A prior based on the luminosity function that disfavors the likelihood of obtaining 
a high redshift for high observed flux is used, such that the function to be minimized according 
to the redshift and the SED type T is: 




-2log(P(M B \T,z)), 



where F \, s is the observed flux, F exp is the expected flux, M a normalization factor, and 
P(Mb\T, z) is the prior based on the luminosity function. 



2.5 Power spectra and dark energy parameter reconstruction 

A fiducial cosmology is assumed and the comoving distance to each galaxy is computed, where 
the same values of the cosmological parameters as the ones used to generate the observed galaxy 
catalog are taken. We lay a grid over the simulated galaxy survey and count the number of 
galaxies in each grid cell n g . The galaxy over-density distribution in cells is then: 

*,(r) = 

n 

where n is the mean number density of galaxies, calculated from the total number of galaxies 
in the simulation divided by the simulation volume. The gridded over-density is then Fourier 
transformed, and the power spectrum is computed as a function of redshift slice. The comoving 
BAO scale is computed using the wiggles- only method. The power spectrum is divided by a 
smooth, wiggle-free power spectrum, and the resulting function is well approximated by an 
empirical slowly decaying sinusoidal function with parameters k B AO, A the amplitude, and a 
decay parameter a: 
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where a is set to 1.4 (see Eisenstein et al^). A chi square minimization is done by varying the 
parameter A and k BA o : 
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where Pg m (k) is the measured power spectrum divided by the wiggle- free reference spectrum 
P re f and Pg m (k, k BA o, A) is given by Eq[lJ The uncertainty on k BA o is given by the solving 
X 2 {kBAO i o~(k,BAo)i A ± o~(A)) — Xmin = 2, 3 at 68% confidence level, where the hat stands for 
the values of the parameters that minimize the chi square. The angular BAO scale 9 S is related 
to the comoving BAO wavelength and the angular diameter distance D A (w, z) at a redshift slice 
z : 

Os = , , (2) 

k B AoP>A(w,z) 

The CMB allows for a calibration of the BAO scale and gives 6 S , therefore we constrain the value 
of Da(w, z). The angular diameter distance depends on the cosmological parameters, such that: 

^ / x c f z dz' 

D A (w,z) = — (3) 

7 ° Jn m (i + z'f + n A (i + z') 3(1+w) 
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Figure 1: Left-hand panel: Reconstructed power spectrum for z = 0.7. Right-hand panel: Evolution of the 
relative uncertainties of the angular diameter distance Da as a function of w and the redshift slice. 

The relative uncertainty on ksAO is equal to the relative uncertainty on Da(w,z) which gives 
the uncertainty on w at redshift z as shown in Fig [I] on the right-hand panel. 

3 Results 

For a very preliminarysimulation of 5 Gpc 3 and z G [0.7; 1.4], we obtain a galaxy catalog of 
2.6 x 10 7 galaxies. The BAO wavelength uncertainties computed by fitting the power spectrum 
shown Fig [l] on the right-hand panel is equal to 3% at z = 0, 7. Looking at the evolution of the 
relative uncertainty on Da as a function of w and the redshift, gives a relative uncertainty on 
the dark energy parameter Aw/w = 13% for ADa/Da = 3%. 
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